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■ Abstract 
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Properties of the horizon mass of hairy black holes are discussed with em- 
^ I phasis on certain subtle and initially unexpected features. A key property 

suggests that hairy black holes may be regarded as 'bound states' of ordinary 

■ black holes (without hair) and colored solitons. This model is then used to 
predict the qualitative behavior of the horizon properties of hairy black holes, 
to provide a physical 'explanation' of their instability and to put qualitative 

I constraints on the end point configurations that result from this instability. 

"q ' The available numerical calculations support these predictions. Furthermore, 

O^' the physical arguments are robust and should be applicable also in more com- 

plicated situations where detailed numerical work is yet to be carried out. 
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I. INTRODUCTION 



In the Einstein-Maxwell-Klein-Gordon theory, static black holes are remarkably simple 
objects. They are completely characterized by their mass and (electric and/or magnetic) 
charge, evaluated at infinity. (For reviews, see As is well-known, this simplicity 

does not persist once non-Abelian gauge fields are incorporated For example, for 

a sufficiently large value of the ADM mass Madm, the Einstein- Yang-Mills theory admits 
distinct spherical, static black hole solutions, labelled by an integer n, all with zero charge 
at infinity^ (see Fig. [^). Although in terms of conserved quantities at infinity these solutions 
are indistinguishable from the Schwarzschild black hole of mass Madm, in the interior their 
geometries and Yang-Mills fields are quite different from one another. Even the horizon 
structure varies as we pass from one n to another; for instance, for ffxed Madm, the horizon 
area is a decreasing function of n. Thus, the rigidity tying conserved quantities at infinity 
to the structure in the interior does not extend to the non-Abelian context. Following the 
common terminology, if this failure occurs we will say that the theory admits 'hair' and refer 
to these solutions as 'hairy black holes'. 

Of central importance to this paper is the notion of the horizon (or black hole) mass 
Mhor in general (Einstein-matter fields) theories with hair. In the Einstein-Maxwell theory, 
for static black hole solutions one can set Mhor = Madm- However, outside these globally 
time-independent situations, this identification is not viable. For, there may be gravitational 
radiation in a region far removed from the black hole, and hence irrelevant for any reasonable 
notion of the horizon mass Mhor, which would nonetheless obviously contribute to Madm- 
Since hairy black holes are known to be unstable, such dynamical situations arise naturally 
in their study. Furthermore, as discussed in Section |1II A| , even in the static context, Madm 
is a poor measure of the mass of these black holes. Therefore, it is desirable to have a 
general, systematic approach to the problem of defining horizon mass. In this paper, we 
will reach this goal by using the recent isolated horizon framework, based on Hamiltonian 
methods @-|Tl[. 

Closely related to hairy black holes are the solitonic solutions to Einstein- Yang-Mills 
equations [0,^,0]. Indeed, it was the Bartnik-McKinnon |]T^ discovery of these solitons 
that shook the then prevailing view that there is no qualitative difference between Einstein- 
Maxwell and Einstein- Yang-Mills theories and paved the way for discovery of hairy black 
holes. The purpose of this paper is to use the isolated horizon framework to gain new, 
qualitative insights in to the relation between colored black holes and their solitonic analogs 
and to point out certain subtleties that arise as we move away from the Einstein-Maxwell 
context and admit increasing complexity in the hair. A key property of the horizon mass 
suggests that a colored black hole can be regarded as a 'bound state' of an ordinary, uncolored 
black hole and a colored soliton. This physical picture turns out to be very useful. It 



^An infinite family of static spherical solutions was shown to exist in [^Q- Furthermore, unlike in 
the Einstein-Maxwell theory, there also exist non-spherical static black hole solutions. Although 
one does not yet have a full control on the entire static sector, the known family of static solutions 
with zero charge at infinity can labelled by the horizon area and two integers Q. 
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ADM Mass 




FIG. 1. The ADM mass as a function of the horizon radius Ra of static spherically symmetric 
solutions to the Einstein- Yang- Mills system (in units provided by the Yang-Mills coupling constant). 
Numerical plots for the colorless, n = 0, and the first two families, n = 1, 2 of colored black holes 
are shown. (Note that the y-axis begins at Af = 0.7 rather than M = 0.) 

enables one to regard the total (ADM) mass as a sum of three terms: the mass of a bare 
(i.e. uncolored) black hole, the mass of the appropriate soliton and the gravitational binding 
energy between the two. Physically expected properties of the gravitational binding energy 
(such as its negative sign) then determine the qualitative behavior of the horizon attributes 
of colored black holes (such as surface gravity). These predictions of the model have already 
been borne out by numerical calculations in simple cases. We will also use the model to 
argue that instability of colored black holes 'stems' from that of solitons and constrain 
the final outcome of small perturbations. Finally, as is well-known, quasi-local definitions 
of mass in general relativity tend not to possess all the properties one is accustomed to 
from field theories in flat space. We will see that while the isolated horizon framework 
provides an internally consistent viewpoint and strategy, this very consistency implies that 
the resulting horizon mass has certain counter-intuitive features. They turn out to have 
important physical implications, in particular to the stability properties of magnetically 
charged Reissner-Nordstrom solutions. 

The transition from Abelian to non-Abelian gauge fields is accompanied by the emergence 
of new scales with dimensions of energy (or length). The classical Einstein-Maxwell theory 
has no such scale. The coupling constant g of the Einstein- Yang-Mills theory provides an 
energy scale Mq = l/(g^/G) (in the c=l units) which signals the onset of new branches 
(labelled by n) of spherical static solutions (see Fig. |I|). In the Einstein- Yang-Mills-Higgs 
theory, the mass and the self-coupling constant of the Higgs field provide additional scales. 
Typically, as the number of such scales increases, the departure from the Einstein-Maxwell 
behavior becomes more significant and the number of unexpected, subtle features grows. 

In Section |T| we recall the definition of the horizon mass in the Einstein-Maxwell theory. 
In Section |T| we consider Einstein- Yang-Mills black holes, discuss a key relation between 
the horizon mass of static, hairy black holes and the mass of the corresponding solitons, 
introduce the heuristic model of colored black holes mentioned above and derive some of 
its predictions. In particular, we will show that the notion of the horizon mass can be 
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used to 'explain' why the magnetically charged Reissner-Nordstrom solutions are unstable 
in the Einstein- Yang-Mills theory [14| although they are stable in the Einstein-Maxwell 
theory For definiteness, in this section we will restrict ourselves to Einstein- Yang- 

Mills theory. However, modulo certain subtleties, most considerations of this section will 
apply also to hairy black holes in more general theories. In Section ^ we consider such 
theories and discuss the relevant subtleties. These arise from the fact that, while in the 
Einstein- Yang-Mills theory hairy black holes belong to disjoint families carrying discrete 
labels (e.g., integers), in more general theories these families intersect. Section provides 
a brief summary. 



II. ISOLATED HORIZONS AND HORIZON MASS 
IN EINSTEIN-MAXWELL THEORY 

To define the horizon mass, it suffices to consider the so-called 'weakly isolated horizons'. 
The precise definition and a detailed discussion can be found in [llj but will not be needed 



here. In essence, a weakly isolated horizon A is an expansion-free, null, 3-dimensional 
sub- manifold A of the 4-dimensional space-time {M,gab), equipped with an equivalence 
class of future-directed null normals which Lie-drag (certain components of) the intrinsic 
connection on A.0 The notion extracts from the definition of a Killing horizon just those 
properties local to the horizon that are needed to obtain a consistent Hamiltonian framework, 
define the horizon energy, angular momentum and other charges, and prove the zeroth and 
first laws [nlJl^. In the non- rotating case now under consideration, the vector field 



on A plays the role of the static Killing field. However, as explicit examples and general 
constructions show, the space-time need not admit a Killing field in any neighborhood of A 



To define the horizon mass Ma, one proceeds as follows. Recall that, in the standard 
asymptotically fiat context, the ADM energy -^adm arises as the boundary term at spatial 
infinity in the expression of the Hamiltonian generating evolution along an asymptotic time- 
translation t"". Now the idea is to use as phase space the sector of general relativity consisting 
of space-times which are asymptotically fiat at spatial infinity and admit an isolated horizon 
A as internal boundary. The Hamiltonian generating an appropriate time-translation then 
has two surface terms, one at infinity, , and one at the internal horizon boundary, E^^ . 
As one might expect, the first is precisely the ADM energy E^^ = -Eadm- The second, E^^ , 
is interpreted as the horizon energy. 

However, in the detailed implementation of this strategy a subtlety arises which will play 
an important role in this paper. At infinity, all metrics in the sector under consideration 
approach a fixed Minkowskian metric rjab- Therefore, we can simply fix a time-translation 



^Here, i"" is equivalent to i"" if an only if = where c is a positive constant on A. Because the 
Maxwell stress-energy condition has the property that —Tah^^ is future-directed, the Raychaudhuri 
equation equation implies that the requirement that (Aj^"^) be expansion-free is equivalent to the 
requirement that i"" Lie-drags the intrinsic (degenerate) metric on A. 
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Killing field of rjab and demand that agree with that vector field in a neighborhood of 
infinity. At the internal boundary, by contrast, the geometry is not universal. Therefore, 
given a candidate time-translation near A for one space-time in our sector, it is not a priori 
clear how to choose the same time translation in another space-time (also belonging to the 
sector). To compute the Hamiltonian, one must fix this ambiguity, i.e., choose appropriately 
a vector field for each space-time in the phase space. In the non-rotating case, it is natural 
to set = ci°' on A. However, while c is constant on A in any one space-time, it may vary 
from one space-time to another, i.e., can be a function (of the horizon parameters) on the 
full phase space. Put differently, in the numerical relativity terminology, one is led to allow 
live evolution vector fields t°- (or, equivalently, lapse-shift pairs that depend on the point on 
phase space) on space-time. 

Let us fix such a live vector field t". The evolution along it defines a vector field Xt on 
the whole phase space. It is then natural to ask if this evolution is Hamiltonian, i.e., if it 
preserves the symplectic structure. Somewhat surprisingly, the answer is in the affirmative 
if and only if there exists a function E^^ of the horizon area oa and electric charge Qa such 
that 

=^^«A + $(i)5gA. (1) 
OTTCt 



for arbitrary tangent vectors 6 to the phase space Here, = stands for equality restricted 
to A, the surface gravity is defined by via f^Vat^ — i^{t)t^ ■, and is the electric 
potential defined by on A. {K{t) is called 'surface gravity' because = ci"', with c a 
constant on A. The potential is unimportant for the purposes of this paper.) Note that 
(0) implies that and are also functions only of oa and Q^; they can not depend on 
other properties of the solution under consideration. 

If is interpreted as the horizon energy, (|I]) is precisely the first law of black hole 
mechanics! Thus, the necessary and sufficient condition for the evolution generated by a 
live vector field t° to be Hamiltonian is precisely that the first law hold. A live vector field 
satisfying this condition is said to be permissible. The Hamiltonian H^^^ generating evolution 
along a permissible vector field t° is given by 

F^*) = f (constraints) d^x + E^H^^ - E^^ (2) 

where S is any partial Cauchy slice extending from the isolated horizon A to spatial infinity. 

One can now ask if such live vector fields exist on the whole phase space .0 The answer 
is in the affirmative [TT|]. Choose any (regular) function of a a and Qa- Then, one can 
show that there exists a permissible live vector field with H(^t) = i^o- Thus, in fact, there 
exists an infinite family of permissible live vector fields, each defining a horizon energy E^ 
and leading to a first law. While the phase space meaning E^^ is transparent, for a general 



^Since tfie volume term in the expression (^) is a linear combination of constraints, only the 
boundary values of the live evolution vector field matter on 'on shell'. Therefore, while discussing 
evolution vector fields, we will be concerned only with their boundary values. 
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permissible vector field t", its space-time significance is not so clear. Is there perhaps a 
canonical choice for the boundary value of t"' on A? In the Einstein-Maxwell theory, the 
answer is in the affirmative. Let us demand that the required live vector field be such 
that in each static (i.e. Reissner-Nordstrom) solution coincides with that static Killing 
field which is unit at infinity. This is a natural choice. Furthermore, it fixes K{to) uniquely 
as a function of oa and Qa- Therefore, although we fixed the normalization of only on 
the static solution, since there exists a unique static black hole solution for each (oa, Qa), 
the functional form of K(^to) is completely fixed on the whole phase space. This in turn selects 
the boundary value of on A and determines everywhere on the phase space. 
It is then natural to set the black hole mass 

Mho. = 4*°^ (3) 

for any space-time in the phase space. This definition of Mhor has several pleasing features. 
First, because to is permissible, (in an obvious notation) the corresponding first law (|1]) takes 
the familiar form 

5Mhor = SaA + ^6Qa , (4) 

but now on the entire phase space, including non-static space-times admitting gravitational 
and electromagnetic radiation. Second, if the evolution vector field coincides with a static 
Killing field, general phase space considerations imply that the total Hamiltonian is constant 
on each connected component of the static sector. Since in the Einstein-Maxwell theory, 
there is no natural energy scale, this constant is necessarily zero. Eq (|]) now implies that on 
static solutions, Mhor = Madm- Finally, in presence of radiation, assuming that the horizon 
extends all the way to future time-like infinity and that the structure of is 'the same 
as that in Reissner-Nordstorm space-times', one can show 

Madm - Mho. = Ej+ , (5) 

the total energy radiated across future null infinity, X+. (Equivalently, even in presence 
of radiation, Mhor is the future limit of the Bondi energy as one approaches along X"*".) 
Thus, in the Einstein-Maxwell theory, the horizon mass Mhor is defined everywhere on the 
phase space and has the properties one intuitively expects. 

Let us summarize. At infinity, Madm is the norm of the asymptotic 4-momentum. 
Because the horizon is in a strong field region, we do not have a 4-dimensional translation 
group on A, whence there is no well-defined notion of horizon 4-momentum. Therefore, 
we need an alternate extra structure to define Mhor- The idea is to use a preferred time- 
translation and identify Mhor with the horizon energy E^^° defined by this translation. 
This task could be carried out because, thanks to the no-hair theorems, we could choose 
a preferred, live evolution field for '^^^ space-times in the full phase space and use it to 
introduce a canonical definition of Mhor on the whole phase space. 



III. YANG-MILLS HAIR, HORIZON MASS AND STABILITY 

This section is divided in to four parts. In the first, we recall the known results from 
|r^,|TT|. In the second, we introduce a model of colored black holes as bound states of the 
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ordinary black holes and appropriate solitons and show that it predicts all the qualitative 
features of the horizon structure of colored black holes in equilibrium. In the third part, 
we use the model to discuss instability of colored black holes and put constraints on the 
final configurations resulting from small perturbations. In the fourth part, we focus on 
the 'embedded Abelian black holes' and show that the horizon mass of these black holes 
is higher if they are regarded as Yang-Mills black holes than it is if they are regarded as 
Einstein-Maxwell black holes. Thus, because the horizon-mass is a phase-space rather than 
space-time notion, in this case it turns out to be 'theory dependent'. This is a striking 
and, at a first glance, surprising property. It has an important physical consequence, briefly 
noted in [^: using considerations of the first three sub-sections, we will show in detail that 
it 'explains' why these solutions are unstable in the non-Abelian sector although they are 
stable in the Einstein-Maxwell theory. This is another striking example of the usefulness of 
the notion of horizon mass. 



A. Horizon mass of colored black holes 

Let us now consider the Einstein- Yang-Mills theory with gauge group SU{2). In this 
case the Reissner- Nordstrom solutions provide a 2-parameter family of static, spherically 
symmetric, Abelian black holes, labelled by (Madm,Q), where Q is now associated with 
a fixed U{1) subgroup of SU{2). In the Maxwell case, if the potential Aa is a connection 
1-form on a trivial bundle, the magnetic charge vanishes identically (and, on non-trivial 
bundles, it is quantized). In the Yang-Mills case, by contrast, there is a 1-parameter family 
of magnetically charged, static solutions (with a fixed magnetic charge Pq) even on the trivial 
SU{2) bundle [0. These are called 'embedded Abelian solutions' because they are isometric 
to a family of magnetically charged Reissner-Nordstrom solutions and the isometry maps the 
Maxwell field strength to the Yang-Mills field strength. The only difference is in the form of 



the connection. We will return to these black holes in the subsection [III D| . Of direct interest 
to this subsection are the families of 'genuinely' non-Abelian, static solutions. Each family 
is parametrized by one continuous parameter (oa)- An infinite number of families, labelled 
by two integers, is known to exist of which spherically symmetric solutions constitute a 
sub-class labeled by one integer of Fig. [0,0. 

We will now show that, in the Einstein- Yang-Mills theory, the ADM mass fails to be a 
good measure of the black hole mass even in the static sector. Consider, for instance, the 
branch of spherical, static black holes labelled by n 7^ (Fig. |ip. Let us decrease the horizon 
area along this branch. In the zero area limit, the solution is known to converge point- wise 
to a regular, static, spherical solution, representing a Einstein- Yang-Mills soliton |]T^,|TB . 
This solution has, of course, a non-zero ADM mass M^^-^, which equals the limiting value 
of M^Qj^. However, in this limit, there is no black hole at all\ Hence, this limiting value 
of the ADM mass can not be meaningfully identified with any horizon mass. By continuity, 
then, M^Qf^ can not be taken as an accurate measure of the horizon mass for any black hole 
along any n 7^ branch. 

Can one use the isolated horizon framework to define Mhor of static black holes in the 
n 7^ families? The answer is in the affirmative but there is a subtlety. Considerations of 
Section |I| leading to eqs (|l|) and (0) go through also in this case [|lT| . Permissible vector 
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fields can be constructed as in Section O and, for eacli permissible vector field we again 

— (t) 

have a well-defined notion of the horizon energy E)^ on the entire phase space, and the 
associated first law. However, since the no hair theorem now fails, we can no longer single 
out a preferred live vector field on the full phase space by using static solutions and 
set Mhor = -Ehor- We could focus just on the electrically charged, static Abelian solutions. 
Surface gravity k^^"^^ of the properly normalized Killing field of these solutions does provide 
us, as in the Maxwell case, with a specific function fi;^^''^(aA, Qa) and we can again use it 
to obtain a preferred, permissible, live vector field tAbei ^^"^ corresponding horizon energy 
^Abei_ Pqj, ^]2ese static solutions, we can again set Mhor = E^"^^. However, n^^"^^ fails to agree 
with the surface gravity of the properly normalized Killing field in hairy static solutions. 
Hence there is no reason to interpret E^^^ as the horizon mass of the hairy solutions. To 
summarize, because the uniqueness theorem fails — i.e., because of the hair — we can not 
single out a canonical function kq of the horizon parameters which agrees with the surface 
gravity of the (properly normalized) Killing field on all static solutions. Therefore, the 
Einstein-Maxwell prescription for defining the horizon mass Mhor on the full phase space 
does not extend to theories with hair. 

However, there does exist a natural strategy to assign a horizon mass to any static 
solution [^. For definiteness, let us consider the 1-parameter family of static, spherical 



solutions labelled by uq. (However, our discussion and conclusions apply to other static, 
hairy black holes as well.) Since they have zero electric charge, along this branch the surface 
gravity (defined by the properly normalized Killing field) is a function /t(no)(oA) only of the 
area a a- By setting hq = K(^no) iii the construction outlined in the last sub-section, we obtain 
a permissible, live vector field f^^ on the entire phase space. The corresponding horizon 



energy E^°'^ is given by [p!0| , p!T 



4"°^ = ^/^''"Ano)Mdr, (6) 

where i?A is the horizon radius (i.e., the horizon area is given by oa = AitR^), and 
/3(„o)(i?A) = (2-Ra) {i^{no){RA))- While this horizon energy E^^°^ has no obvious space- 
time interpretation for general space-times, for static solutions on the branch n = no, it is 
natural to set 

1 r^A 



(i?A) = E^-^\Ra) = ^ // Ano)(r)dr . (7) 



We will do so. Note that M]^"°'*(_Ra) vanishes at -Ra = (and then monotonically increases 
with area). Hence this definition of the horizon mass is free of the problem faced by the 
ADM mass, discussed above. 

Finally, for any n one can relate the horizon mass m£J to the ADM mass of static black 
holes. Recall first that general Hamiltonian considerations imply that the total Hamiltonian 
(§) is constant on every connected component of static solutions (provided the evolution 
vector field agrees with the static Killing field everywhere on this connected component) 
PJTT[|. In the Einstein-Maxwell case, there is only one connected component and there is no 



energy scale in the theory. Therefore, we could conclude if(*o) = on all static solutions. 
In the Einstein- Yang-Mills case, by contrast, there is an infinite number of connected com- 
ponents of static solutions and there is a natural energy scale, Mq = l/gy/G. Therefore, 
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on each distinct connected component, in principle, the Hamiltonian could be a different 
multiple of Mq. This is exactly what happens. Since the Hamiltonian is constant on any 
n-branch, we can evaluate it at the solution with zero horizon area. This is just the soli- 
ton, for which the horizon area oa, and the horizon mass M^or vanish. Hence (0) implies 
ff{to,n) ^ ^in)_ rj.^^^^ conclude: 

M^:^ = mSm - mS (8) 



on the entire nth branch [|TO|,|rT|. Thus, the ADM mass contains two contributions, one 
attributed to the black hole horizon and the other to the outside 'hair', captured by the 
'solitonic residue'. 



B. A physical model of colored black holes 



Considerations of Section |III A| lead to a heuristic, but quite powerful physical model of 
a colored black hole. In this subsection we will develop this model and derive some of its 
consequences on the equilibrium properties of static black holes. 

Recall that for static black holes in the Einstein- Maxwell theory, the future limit M^^^^^ 
of the Bondi mass along X"*" equals the horizon mass M^or- What is the situation for the 
non-Abelian black holes of the Einstein- Yang-Mills theory? In the zero area limit along 
any n 7^ branch, one is left with a soliton, i.e., with a regular bound state which persists 
all the way to and thus contributes to the mass at Because the solution is static, 
the ADM mass, the Bondi mass and the mass at are all equal. As we move away from 
zero area, a true black hole appears. For small areas, implies that M^°'' is small. But, 
in relation, M^^}^ can be large because 'there is still a solitonic contribution to the ADM 
mass', quantified in (|]). In the Einstein-Maxwell case, there are no solitons; there is no 
'residue' at apart from the black hole itself. This is why for static black holes, Mhor 
equals Madm and charges at infinity suffice to characterize the fields on the horizon. In the 
Einstein- Yang-Mills case, by contrast, the conserved quantities at infinity do not in general 
capture the whole story. For hairy black holes, the structure at X"*" (or i^) does not rigidly 
determine the structure at the horizon; there is a solitonic residue at which accounts for 
the lack of rigidity. This suggests that we think of a colored black hole as a 'bound state ' of 
an ordinary, 'bare' black hole and a 'solitonic residue'. Indeed, we can trivially rewrite the 
total mass of (|^) as a sum of three parts, one referring to the soliton, one to the 'bare' black 
hole and one to the binding energy: 

= + M^°; (i?A) + ^binding (Ra) , (9) 

where M^^{Ra) = Ra/2G is the horizon mass of the Schwarzschild black hole of radius Ra- 
Thus, one can say that the effect of the solitonic residue is to 'dress' the black hole mass 
from its 'bare value' M^^{Ra) to M]^"](i?A)- The fact that the total space-time mass Madm 
is just a sum Madm = +M^^i , without an explicit interaction term, can be understood 
as follows. Recall that, in general relativity, the total Lagrangian or the Hamiltonian is just 
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FIG. 2. The horizon mass Mhor as a function of the radius Ra (in units provided by the 
Yang-Mihs coupHng constant). Numerical plots for n = 0, 1, 2 spherical, colored black holes in the 
Einstein- Yang-Mills theory are shown. 

a sum of the gravitational part without matter and the matter contribution. There is no 
explicit interaction term between matter and gravity; the interaction is 'absorbed' in the 
description by simply replacing the flat metric in the 'bare' matter terms by the curved, 
physical metric, making the full expression covariant. In an analogous way, the interaction 
between the soliton and the black hole is expressed covariantly by simply replacing the 'bare' 
black hole mass Mj^^lR^) by the 'dressed' mass M^l{R/^) as in (H). 

Let us adopt this model and work out some of its implications. More precisely, the idea 
is to use this model to derive quahtative features of the properties of colored black holes 
assuming, whenever necessary, properties of colored solitons. Since the colored black holes 
solutions have been obtained numerically, we do not have simple analytical expressions of 
their horizon attributes such as surface gravity and horizon mass. Furthermore, typically 
such quantities have been computed only for low values of n. In the remainder of this 
sub-section we will show that all the qualitative features of the known numerical plots of 
these quantities as functions of n and horizon radius can be derived from our model. 
Furthermore our predictions apply to all values of n and -Ra- 

First, since -^binding is the gravitational binding energy between the soliton residue and 
the bare black hole, it must be negative. Thus, the model together with the elementary 
physical fact about the gravitational binding energy leads us to the first prediction: 

i) M]^Q^(i?A) < M^2{R/^) for all n > and all values of the horizon radius R/^. 

At first this conclusion seems somewhat counter-intuitive since, for a given horizon area, one 
normally thinks of the colored black holes as 'excited states' and the n = 0, Schwarzschild 
black hole as the 'ground state'. However, that intuition comes from the properties of the 
ADM mass: it is known that M^Qj^(i?A) > ^adm(-^a)- Explicit numerical calculations 
have in fact confirmed i) for low values of n (see Fig. 0): while M^j^f^{R/^) does grow with 
n, M^^i grows so that i) holds. The model predicts that this inequality will hold for all 
n > 0. 



10 




o 







2 



4 



6 



FIG. 3. Numerical plots of = 2k(„) i?A as function of the horizon radius Ra for n = 1,2 
spherical, colored black holes in the Einstein- Yang-Mills theory. For large Ra, all curves asymp- 
totically approach the /3 = 1 line. 



for all n and Ra- Hence, it follows that [P(n) — P{o)) (Ra) < 0, whence we conclude 



Again, this prediction is borne out by detailed numerical simulations for low n (see Fig. 

We can now use the physically expected properties of the gravitational binding energy 
to make sharper predictions. Let us fix Ra and vary n. Then, the mass M^J^ of the bare 

(n) 

black hole is fixed but, as is well-known, the soliton mass M^^j increases with n. Now, since 
the absolute value of the gravitational binding energy increases as the mass of either of the 
two bound objects grows, | -^binding I must increase with n. This implies: 

iii) For any fixed value of Ra both the horizon mass M^^l and the surface gravity K{n) are 
monotonically decreasing functions of n. 

As Fig H and Fig ^ show, these predictions are borne out numerically for low n. 

Next, let us consider the complementary situation and keep n > fixed and vary Ra- 
Now the soliton mass m'^^x is fixed and the bare black hole mass M^^^{Ra) = Ra/'^G in- 
creases with Ra- Hence, again (the absolute value of) the gravitational binding energy must 
increase with Ra- Thus, we conclude: 

iv) For any fixed n, {Ra/2G — mI^I){Ra) is a monotonically increasing function of Ra 
whence, from (Qj, P(n){RA) < 1 for all n > and Ra- 

Again, Fig. ^ and Fig. ^ show that this prediction is borne out by numerical studies for low 
values of n. Furthermore, this property can be directly derived from the definition (^ of 



Next, let us use the expression ([7|) of the horizon mass. From i), we conclude: 




ii) [i^{n) — 1^(0)) (Ra) < for all n > and Ra- 
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m1^2{Ra) and prediction ii) above. Indeed, 

= M£i-M£ = /^''"(K(o)-KH)(r)dr 

and ii) ensures that the integrand is positive. 

Next, we can put a lower bound on We will use the fact that, for any fixed n the 

ADM mass M^^-^IRa) is a monotonically increasing function of Ra- Since M^^^ is fixed on 
the n-th branch of static solutions, it follows from (M) that 



v) mI^I{Ra) is also a monotonically increasing function of Ra- (0j now implies that, for 
all n and Ra, (3{n) o-nd are positive functions of Ra- Hence M^^ is also positive for all 
n and Ra, except for Ra = where it vanishes - 

Finally, the model has a prediction on the asymptotic behavior (for large Ra) of mI^I{Ra), 
P{n){RA) and K(„)(_Ra)- Using the known property of the ADM mass, M^-^IRa) > 
^adm(-^a) = -Ra/SG, we conclude: 



Hence, using prediction i) we conclude that, for any given n, the curve M£(i?A) hes between 
the two parallel lines /i(-Ra) = -Ra/2G — M^^^ and /2(-Ra) = -Ra/2G. Furthermore, since 
mI^I{Ra) is a monotonically increasing function of Ra, it follows that: 

vi) For any n, the curve showing the dependence of M^^^{Ra) on Ra becomes asymptotically 
parallel to the two lines which hound it, i-C-, has slope 1/2G for large Ra- Hence, by the 
curves showing the dependence of I3(n) on Ra asymptotically approach the curve (3{q){Ra) = 1 
and the curves showing the functional dependence of of Ra asymptotically approach the 
curve K(o)(-Ra) = l/2i?A- 

These six properties account for all qualitative features of Fig. | and Fig. || i.e., the 
dependence of M^l and on n and Ra- In particular for Fig. ||, we have shown that our 
model implies that the curves: a) intersect the y-axis at distinct points between and 1 
and lie in the upper right hand quadrant, b) the curves never intersect, c) have the property 
that higher the n, lower the curve, and d) for large Ra become asymptotically tangential to 
the curve P(o){Ra) = 1. To arrive at this qualitative understanding of the horizon properties, 
in addition to the model, we used only two known qualitative properties of the ADM mass 
||T8| : Mgj^j-' increases with n and, for any given n, M^^-^^ increases monotonically with Ra- 



C. Instability of colored black holes 

We will now argue that the model also 'explains' the instability of the colored black 
holes in terms of the instability of the solitons. Under small perturbations, the solitons are 
unstable; the energy 'stored' in the bound state of the soliton is radiated away to future null 



infinity, X"*" |jT9[. In our model, each static black hole with n 7^ is accompanied by a soliton 



which 'hovers around all the way to i"*"' and carries information about its hair. Therefore, it 
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FIG. 4. An initially static colored black hole with horizon Ain is slightly perturbed and decays 
to a Schwarzschild-like isolated horizon Agn , with radiation going out to future null infinity . 

is natural to expect that under small perturbations in the exterior region, this soliton would 
be unstable and radiate away across X+ and A. At the end of the process, there should be 
no solitonic residue at The end product should thus be an uncolored (i.e. Schwarzschild) 
black hole whose horizon mass equals the future limit Mj+ of the Bondi mass along 
X"*". Our suggestion that the instability of the colored black hole is due to the instability of 
the accompanying soliton is supported by the fact that all the colored black holes on 
the nth branch have the same number (namely, 2n) of unstable modes as the nth soliton. 
Of course the detailed features of these unstable modes will differ especially because they 
are subject to different boundary conditions in the two cases. 

In general, however, even if one component of a bound system is unstable, the total 
system may still be stable if the binding energy is sufficiently large. An obvious example is 
provided by the deuteron nucleus — a bound state of a proton and a neutron. Therefore, to 
see if our bound state of a bare black hole and a soliton is in fact unstable, let us examine 
the energetics. If the decay suggested above is to occur, the total energy in the initial 
bound state, ^^i^itiai = ^^(/^^^"^i) + M^^, should exceed the energy in the final black hole 
E'finai = ^hor(-^A"^')- Siucc the area of the final uncolored black hole is at least as large as 
the area of the initial colored black hole, the decay can occur only if the following inequality 
holds: 

Mtl (i?A'"^') + > m£; (b:^'-') . (10) 

The inequality is in fact satisfied. For, the left side is the ADM mass M^^-f^{RA) while 
the right side is the ADM mass M^^yi{R/^) of the uncolored black hole of the same radius 
^mitiai^ and the ADM mass is known to increase with n. Thus, the energetics are such that 
the instability of the soliton can persist also in the bound state, forcing the solitonic residue 
to radiate away. 

Let us therefore assume that the process illustrated in Fig. |^ does take place and examine 
its consequences. Thus, we assume that the event horizon has an initial isolated component 
^initial with color n and radius i?^'*'^' and there is a small perturbation on the initial partial 
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Cauchy surface M. The perturbation evolves and, because of the instabihty of the sohton, 
the sohtonic residue is radiated away and we are left with a final, uncolored isolated horizon 
Afinai with radius B^^^ (satisfying R^^^ > i?^'*'^') and energy Ej+ radiated across future 
null infinity X"'".^ The ADM mass of the space-time does not change in this process. Hence 
neglecting the energy in the initial perturbation (i.e., in the limit 'from above' of a sequence 
of perturbations with decreasing energy) we have: 

= M£(i?r^^) + Mi") = M£;(i?r ) + . (n) 

Using (^ and the expression of the ADM mass of the uncolored (i.e. Schwarzschild) black 
hole in terms of its horizon radius, this equality can be re-expressed as 



I TTiinitial 

sol I binding 



^ - + . (12) 



2G 

Thus, the 'available energy' for the process is given by: 

_ A/f(")_ I pinitial I (-[rCS 
-^avail ~ ^'^sol I -'^binding I • ) 

Note that -E^vaii "^^^ t)e computed knowing just the initial configuration. However, the distri- 
bution of this energy in to a part which goes in to increasing the horizon area and the part 
which gets radiated to null infinity will depend on the details of the initial perturbation. 

Our model enables us to make further qualitative predictions about this process. Note 
first that, since for a fixed n, the binding energy | -EbJnifng I is a monotonically increasing 
function of i?^'*"^^ and M^^^ is independent of -Ra\ the available energy -Eavaii decreases as 
^mitiai increases. In this sense, for a given n, a larger colored black hole is less unstable than 
a smaller one. This prediction is supported by the fact that for the n = 1 colored black 
holes, the frequency of all unstable modes is a decreasing function of the area, whence the 
characteristic decay time grows with area [0,^. Next, let us fix i?^'*"^^ and let n vary. 

Then, since the available energy can be re-expressed as -Eavaii = ^adm ~ ^MiMi it follows 
that it increases with n. In this sense, for a given i?^'*^^', higher the n, more unstable the 
colored black hole. This prediction is supported by the fact that the number of unstable 
modes is given by 2n. It would be interesting to analyze if there exist further, more detailed 
correlations. 

Finally, it is natural to ask if the horizon radius necessarily grows in this process. Using 
an additional input, we will now show that the answer is in the affirmative. Recall from 
IO| , p!T[] that each isolated horizon in the Einstein- Yang-Mills theory carries a well-defined 



magnetic charge Pa^ 



Pa = -i- ^ I F I e (14) 



^In a realistic situation, in distant future the event horizon will become isolated only asymptot- 
ically and will probably not admit a finite isolated piece Agnai- However, our conclusions will 
continue to hold in this case if the approach to isolation is sufficiently fast. 
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where the integral is taken over any 2-sphere cross-section of A, e is the area 2-form on 

1 

this cross-section, and | F \= \Ti{Fab^°'^) (-Fc^e'^'^)] 2 . The additional fact we need is that Pa 
vanishes if and only if n = 0. Therefore, in the process under consideration. Pa necessarily 
changes. Let us now assume P^'*^^' = R^'^^ and arrive at a contradiction. Since the 
expansion of the null normal £° of any event horizon is necessarily non-negative, it follows 
that the expansion is necessarily zero. The Raychaudhuri equation then implies Tahd-^'v^ = 
for all vectors f " tangential to the event horizon and the form of the Yang-Mills stress-energy 
Tab then implies that the field strength must satisfy Fab^°'v^ = 0. One can now use the 
reasoning of |]1T| (see Eq (VI. 7) to Eq(VI.9) of that paper) to conclude that Pa is conserved 



all along the event horizon (even though some of the isolated horizon boundary conditions 
may be violated in the intermediate region). This is however impossible because while p^'t^'^i 
is non-zero P^^^ vanishes because the end product of the process is a Schwarzschild black 
hole. Hence our assumption P^^*'*^^ = R^^^ can not hold; the area of the event horizon must 
necessarily increase in this process. Hence (|1^ implies: 

Ei = (il^o"^- I E^^^ 1) - ^ (Pl-i - Pr"^i) < . (15) 

There may well exist a more sophisticated versions of this argument which provide stronger 
bounds on the change in the horizon area in terms of P™t'^i and P^^*''^^ and hence a better 
upper bound on the energy that can be radiated away to X+. 



D. Instability of Embedded Abelian solutions 

Let us now consider the 'embedded Abelian solutions' with magnetic charge. In the 
Einstein-Maxwell theory, the electromagnetic connection is defined on a non-trivial U{1) 
bundle of a fixed Chern-class, corresponding to the fixed value Pq of the magnetic charge 
while in the non-Abelian case, the Yang-Mills potential is defined on a trivial SU{2) bundle. 
(In the non-Abelian case, Pq has the same dimensions as the inverse 1/g of Yang-Mills 
coupling constant. In the conventions/units generally used in the literature, Pq = ±1.) In 
either case, a detailed analysis shows that Eqs (1) and (2) continue to hold, without any 
extra terms. Since the presence of an electric charge plays no essential role in this analysis, 
for simplicity, let us focus on the sector of the phase space where it vanishes. Then, in both 
cases, we are left with a 1-parameter family of static solutions which can be coordinatized 
by the value of the horizon area. Using the expression of the surface gravity 

- ^ I' - pTJ 

of this family, we can introduce a permissible evolution field to and obtain the corresponding 
first law: ^ 

Let us integrate this law to obtain the horizon Energy: 

4'°'(flA) = ^(l + ^)H-C (16) 
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where C is the integration constant. For colored black holes, the constant was determined 
1^,11] by requiring that E^^°^ should vanish in the limit R/\ goes to zero. For (embedded) 
Abelian black holes under consideration, by contrast, the minimum value of i?A is VG \ Pq \ 



when the black hole becomes extremal. Therefore, we need a new strategy to determine the 
constant C. 

In the Einstein-Maxwell case, it is natural to appeal to the duality invariance of the 
theory. The sector of the theory with zero magnetic charge but arbitrary electric charge Q, 
we know P,pl]| that the horizon energy is given byQ: 



Hence, by duality invariance, we conclude C = 0. 

One could also have reached the same conclusion by starting ab-initio in the dual descrip- 
tion, i.e., by considering from the beginning the vector potential for the electric field as the 
configuration variable and the magnetic field as the momentum variable and allowing the 
magnetic charge P to be arbitrary but setting the electric charge to zero from the beginning. 
One could then work on a trivial bundle and, using the dual of the standard framework of 



Tl| , conclude that C = 0. To summarize, because of our choice of the evolution vector field 
tg, we can interpret as the horizon mass M™ of the Einstein-Maxwell magnetically 
charged black holes and each of the two arguments given above leads us to the result 



Let us now turn to Yang-Mills theory. To determine the constant C in we can 

no longer appeal to duality invariance or use a dual description with the vector potential 
of the electric field as the configuration variable. Thus, we need a new input. Let us first 
consider the colored, spherically symmetric, non-Abelian black hole space-times discussed 
in the last sub-section. It is known that as n tends to zero, these colored, non-Abelian black 
hole space-times tend to the magnetically charged, embedded Abelian black hole space-time 
provided the horizon radius satisfies the bound Ra > VGPo If Ra = VG | -Pq I, the 



magnetically charged black hole is extreme. The new input is that as n tends to infinity, 
/3(n){RA) tends to zero for < -Ra < VGPq (see Fig. |]). Hence, using (^ it follows that as 
n tends to infinity, keeping R^ = VGPq, the horizon mass M£(i?A) tends to zero. Hence, 
in the Yang-Mills phase space, we are led to set the constant C of (^) to: C = —Po/\/G. 
Consequently, in Yang-Mills theory, the horizon mass of the embedded Abelian solutions is 
given by |1T0|: 



^In the reasoning used in [PJTT||, a key role is played by the fact Q is allowed to vary, i.e., that the 
first law reads: (5Ma = (K/87rG)5aA + ^5Q, and by the fact that on this phase space with variable 
Q, there is no constant of dimensions of energy we can construct from the only available constants 
G and c. 
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2g[ ^ Rl I VG 



<or = ^ (19) 



Thus, the horizon mass in the Yang-Mills theory is lower than that in the Einstein theory. 

It is this difference that accounts for the strikingly different stability properties of these 
black holes in the two theories. Since the solutions are asymptotically flat at spatial infinity 
in the standard sense, their ADM mass is determined only by the asymptotic behavior of 
the space-time metric and is given by: 

in both theories. Thus we conclude: 

Madm-M™ = 0, while Madm-M™ = ^ (21) 

V G 



Thus, as guaranteed by general arguments from symplectic geometry, the difference between 
the ADM and the horizon mass is indeed constant on the connected component of magnet- 
ically charged static solutions but the value of the constant is theory dependent. In the 
Einstein-Maxwell theory, the value is zero. All the energy in the solution is associated with 
the horizon and there is no provision in the 'energy budget' for energy to be radiated away. 
In the Einstein- Yang-Mills theory, the energy budget does have such a provision. The future 
limit of the horizon mass (= M^^^) is less than the future limit of the Bondi mass along 
X"*" (= Madm)- Therefore, there is a mismatch at Although in the limit n ^ oo there is 
no regular soliton,^ since 

limM^J ^0, 



from the standpoint of energetics, there is still a solitonic residue at Hence, as a limit- 
ing case of the situation considered in the last sub-section, our reasoning implies that the 
embedded Abelian black holes should be unstable in the Einstein- Yang-Mills theory even 
though their Einstein-Maxwell counterparts are stable. Detailed analysis |T^ has shown 



that this is precisely what happens! Furthermore, we can compute the 'available energy' for 
the black hole of radius Ra in EYM. It has the form E^™ = Pq/2R^. This suggests that 
the frequencies of the unstable modes will be smaller for larger black holes. 

It may seem surprising at first that there is such a marked difference in the stability 
properties in spite of the fact that the space-time geometries of these Einstein-Maxwell 
and Einstein- Yang-Mills black holes are identical. However, since these geometries 'live' in 
distinct phase spaces, modes which can trigger instability within one phase space can be 
absent in the other. Indeed, stability is not a property of a specific solution in isolation but 
refers also to other 'accessible' states and the set of accessible states is different in the two 
phase spaces. The notion of the horizon mass can capture this difference because it is a 
phase space notion rather than a space-time notion. 
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It has been shown that for r > 1 the limiting solution approximates the extremal RN solution 



|23| 
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FIG. 5. The ADM Mass as a function of the horizon radius Ra in theories with a built-in 
non-gravitational length scale. The schematic plot shows crossing of families labelled by n = 1 and 
n = 2 at i?A = B^a"'- 

To summarize, these magnetically charged black holes bring out the fact that the horizon 
mass Mhor of an isolated horizon A can not be determined knowing only the space-time 
geometry unless one knows before hand which phase space this space-time belongs to. At 
first, this 'ambiguity' in the horizon mass seems surprising. However, as we have seen, it 
has a deep physical implication. In the next section, we will see that a further subtlety 
associated with the notion of horizon mass arises as we consider theories with additional 
dimensionful constants and more complicated hair. 

IV. HAIRY BLACK HOLES IN MORE GENERAL THEORIES 



Let us now turn to more general theories with hair, by allowing Higgs or Procca fields in 

These 



addition to the Yang-Mills fields, or considering Einstein-Skyrme theories [24,25,18 



theories have additional dimensionful constants which trigger new phenomena. The static 
black hole sector of such theories enables one to analyze these novel aspects of the interplay 
between gravity and the nonlinear matter fields. Of primary interest to us is the 'crossing 
phenomena' of Fig. ^ where curves in the 'phase diagram' (i.e. plot of the ADM mass versus 
horizon radius) corresponding to the two families of static solutions intersect. With due 



attention to this phenomenon, considerations of Section |T| will extend to hairy black holes 
in these more general theories. In this section we will discuss the relevant subtleties that 
must be taken in to account. 

We will begin with a qualitative argument to elucidate the 'origin' of the crossing phe- 
nomenon. Let us first ignore gravity and consider these theories in Minkowski space-time. 
The coupling constants and other parameters (such as the vacuum values of fields) in these 
theories provide a length scale which determines the 'size' -R^™'^ of the solitonic solution. 
Let us now couple these theories to general relativity. Gravitational effects, being attractive, 
can only reduce this size. Hence, if a black hole is to exhibit any hair at all, its horizon 

(In fact, using more detailed arguments involving the 



radius -Ra must be less than R^^^. 
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pressure of matter fields near the horizon and at infinity, one can give a sharper bound: 
Ra < {2/3)R^l''^ H.) Thus, in these theories which have an in-built length scale already in 
Minkowski space-time, the horizon radius i?A of hairy black holes is bounded above by that 
length scale which is independent of Newton's constant. Note that, Yang-Mills theory in 
Minkowski space has no such built-in length scale whence there is no bound on the horizon 
radius of hairy, static Einstein- Yang-Mills black holes. 

If there is an upper bound to the horizon radius, it is natural to ask what happens 
in the phase diagram to the individual, connected branches corresponding to static black 
holes. Recall first that we can view the space of static black hole solutions as the manifold 



of extrema of the ADM mass functional on phase space |^. In the Einstein- Yang-Mills 
theory, each connected component of this manifold is a sub-manifold without boundary. Let 
us suppose that the topology of these sub-manifolds does not change abruptly as the new 
coupling (and other) constants are slowly switched on to arrive at more general theories. 
Then, if the horizon radius of the new theories is bounded from above, the sub-manifolds 
which are disconnected in the Einstein- Yang-Mills theory must merge smoothly in these 
theories.^] Then, in the projection of this smooth sub-manifold to the (-Ra,^adm) plane, 
one would see two branches which cross at the point corresponding to the maximum value of 
-Ra- This is indeed what is found in numerical investigations of the spherically symmetric, 
static black hole sectors of Einstein- Skyrme, Einstein- Yang-Mills-Proca, and Einstein- Yang- 



Mills- Higgs theories |24-26| and, more recently, in the static axially symmetric black hole 



sector of Einstein- Yang-Mills-Higgs theory [^. Our discussion suggests that this crossing 



phenomenon will occur in any theory that contains a length scale independent of G and 
allows hairy black holes. 

Let us start by analyzing the way two branches merge at the cross-over point. For 
simplicity, let us focus on static black holes that carry no other charge at infinity beyond 
the ADM mass. As long as the energy momentum tensor in the theory satisfies appropriate 
energy conditions, the static black holes can be assumed to admit a Killing horizon with 



a bifurcation two surface |28|. For such black holes, arbitrary perturbations are known to 



obey the first law of black hole dynamics [21 



SMabm = 5aA . (22) 

OTT 

Now consider the phase diagram (i.e., the plot of the ADM mass as a function of the horizon 
radius) for this family and focus on the intersection point between the two branches. If 
the two curves intersect at a finite angle, the first order change in the mass for a given 
infinitesimal change in the horizon area would depend on the branch used to approach the 



■^In the case of the 'embedded Abelian family', by contrast, the sub-manifold of static solutions 
has a boundary corresponding to extremal black holes. However, in this case, the phase space itself 
has a boundary consisting of space-times whose isolated horizons have zero surface gravity and 
the sub-manifold just ends at this boundary. If sub-manifolds of static colored black holes were to 
terminate without merging, these sub-manifolds would have a true boundary in the interior of the 
phase space. 
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intersection point. In view of (]22|), this would contradict the fact that the surface gravity 
has a well defined value at the solution corresponding to the intersection point. Hence, we 
conclude: 

vi)' at any cross-over point, the two branches must have the same tangent vector. 

Explicit calculations have confirmed this behavior in all theories where crossing phenomena 
has been seen to occur. 

Next, let us examine how the crossing phenomenon affects the notion of the horizon mass. 
As in Section [III A|, for all theories under consideration, given any permissible evolution 



vector field we are again led to a horizon energy E^^ and the corresponding first law (p. 
The delicate point again is the selection of a preferred evolution vector field. Because of the 
presence of hair — more precisely, because the value of the horizon area fails to determine 
a static black hole uniquely — we can no longer use the strategy that was successful in the 
Einstein-Maxwell case to define a canonical evolution vector field and use E^l°^ as the 
horizon mass Ma on the full phase space. However, as in the Einstein- Yang-Mills case, we 
can again focus on any one branch, say the nth, of static solutions, use the surface gravity 
^(n)(-RA) to select a permissible evolution vector field t"^-, and construct the horizon energy 

E^^'^ and set (a la (0)) the mass M^^l of any black hole on the nth branch to be: 

mS(^a) = 4"H^a) = ^ Z^""" 0in){r)dr . (23) 

This strategy is tenable except at the cross-over points. At these points, there is a problem: 
Because the function depends on n, the value of the horizon mass is now ambiguous. 
If branches n and n + 1 cross, as in Section |111 B| we have /?(n+i)(-RA) < P{n){RA), whence 



"^hor*"^'' < ^hor- ^aw in Sectiou P^IID| that, since the notion of the horizon mass refers to 
the phase space, the value of the horizon mass of a black hole is not determined simply by 
its space-time geometry. Here we see another, and starker, aspect of the same phenomenon. 



In the magnetically charged black holes of Section |III D| , the same space-time appeared 
in two different phase spaces and the value of its horizon mass depended on which phase 
space one used. In the present case, there is a single phase space. However, we can use 
either the nth family of static solutions and choose a permissible evolution vector field 
using the functional dependence of the surface gravity on the horizon radius Ra to fix 
the normalization, or we can do the same using the (n + l)th family. Since the functional 
dependences of and K(n+i) on the horizon radius are necessarily different, we obtain two 
distinct vector fields and two distinct Hamiltonians which generate time-evolution along 
them. Each of the resulting horizon energies ^i"^ and E^^'^^' is a well-defined function on 
the entire phase space and leads to a first law (0). Although the two energies differ, there is 
complete consistency because the two energies arise from two different evolutions (i.e., two 
different normalizations of the evolution vector field on the horizon). The ambiguity arises 
only when we want to identify these energies with horizon masses and is confined just to 
the cross-over points. As remarked in the Introduction, definitions of quasi-local mass in 
general relativity invariably encounter unexpected features and the ambiguity in the horizon 
mass that arose (in Section [111 L)| and) here is the counter-intuitive feature encountered in 
the present Hamiltonian approach. 
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However, in spite of — or rather, because of — this ambiguity, one can extract useful 
information from this definition of horizon mass. Let the nth and (n+ l)th branch intersect 
at the horizon radius i?^**^"^. Then, using (|^) it follows that: 

mS]m(^a) -M£(i?A) = MiJ and M^^^\r^) - Mi::'\R^) = M^^-^'K (24) 
Hence, using (|^), at the intersection point we obtain: 

M^Ml^tl = + dr = A^r^) + ^ y/ f^inMr) dr (25) 

The second equality yields: 

^n+i) _ = _ ^ (r) dr - — " ir)dr = —f(3{r) dr (26) 

where the last closed counter integral is performed by first moving along the nth branch 
from i?A = to -Ra = R^^^'^^ , then moving back along the n + 1th branch to R/^ = 
and finally sliding down the y axis to the point R^ = on the nth branch (see Fig.^.0 
Note that the counter integral is non-zero because and /3(n+i) are distinct functions 
of i?A, or equivalently, precisely because the horizon mass undergoes a jump at the cross- 
over point. This equation is striking in that it provides a quantitative relation between the 
horizon properties of hairy black holes and masses of solitons, even though the two belong 
to completely different sectors of the theory. The equality was first checked numerically 
in certain static, axi-symmetric hairy black hole solutions in the Einstein- Yang-Mills-Higgs 
theory |2^. Given the generality of our considerations, one expects this relation to hold also 



for more general matter sources which lead to the crossing phenomena, such as Yang-Mills- 
Higgs, Yang-Mills-Proca and Skyrme fields. 

Finally, our discussion of the qualitative behavior of equilibrium properties (Section [ill Bj ) 



and instability of colored Einstein- Yang-Mills black holes (Section [111 Q) will extend also to 



these more general black holes with two provisos: a) due care is taken of the ambiguity in 
Mhor at the cross-over points; and, b) property vi) in Section [1IIB| regarding the asymptotic 



behavior of horizon quantities for large i?A^ replaced by vi)' noted above on the behavior of 
the phase diagrams at the cross-over points. 

Remark: The above discussion can be used also to shed light on the behavior of the ADM 
mass of solitons in the Einstein- Yang-Mills theory, in particular its dependence on n. Let us 
begin with the collection of Einstein- Yang-Mills-Higgs theories parametrized by the vacuum 
expectation value rj of the Higgs field. Now, the static sector of these theories consists of 
extrema of the ADM mass at a fixed area of the internal boundary |^ and, when f] = 0, the 
extrema are configurations where the Higgs field vanishes identically. Thus, we can consider 
the static sector of the Einstein- Yang-Mills theory as the limit rj —>■ of the static sector of 



Since the 'y-axis' corresponds to zero area, points on it do not belong to the black hole phase 
space considered here. Therefore, one should close the contour by moving along the vertical line 
R/\ = e and then take the limit e — > 0. 



21 



the Einstein- Yang-Mills-Higgs theory. Now, in the Einstein- Yang-Mills-Higgs theories, the 
crossing phenomenon occurs and the masses of sohtons are related by Eq. (|26|) for all values 
of rj. Taking the limit when 77 — > we arrive at an identical equation, now relating the ADM 
masses of the Einstein- Yang-Mills solitons. In this limiting case, the branches do not cross 
at a finite value of the horizon radius. However, from the discussion at the beginning of this 
section, one expects that as rj tends to zero, the value of the horizon radius at the crossing 
point tends to infinity. Now, from the Fig. 0, and more generally from prediction vi), we see 
that -/5("+^))(r) as i?A ^ 00. Thus, in the limit, all the (3 curves meet at infinity, 
including the Schwarzschild branch corresponding to /3(o)(-Ra) = 1- Hence, it is plausible 
that the area bounded by the two curves /5(„) and /3(n+i) in the plot of /5 as a function of 
i?A is finite, i.e., that the contour integral on the right side of (p6|) is well-defined also in 
the Einstein- Yang-Mills theory. Thus, we obtain an explicit prediction for the value of the 
ADM mass of the n-th soliton in this theory: 



Again, the equation is striking because it enables one to compute soliton masses from sur- 
face gravity of black holes, now in the Einstein- Yang-Mills theory. To our knowledge, this 
prediction has not been tested. Note, however, that it uses not only the framework presented 
in this paper, but also a continuity assumption on the embedding of the static sector of the 
Einstein- Yang-Mills theory in that of the Einstein- Yang-Mills-Higgs theories. Nonetheless, 
a qualitative consequence of ( P7[ ) is borne out by numerical simulations. The behavior of 
in Fig. ^ suggests that Mg")^^"* — M^^^ approaches zero as n tends to infinity and this 
behavior is supported by the numerical study of low n soliton solutions [l^ll^ , and by rig- 
orous methods in [0]. It would be interesting to test P7| ) in a more direct ways. Finally, 
note that if we combine Eq. p7| ) with Eq. (H) we arrive at an expression for Madm(-Ra), 



The first term is what we have called the 'bare mass', and the second term could be inter- 
preted as the 'mass of the hair', corresponding to the soliton mass plus the binding energy 
as in Eq. (|^). 



Because of the black hole uniqueness theorems in Einstein-Maxwell theory [0,0 it was 
widely believed that black holes would have no hair also in theories more complicated, non- 
Abelian gauge fields. As soon as this expectation was shown to be incorrect 0,^, the subject 
of hairy black holes became a focal point in the mathematical physics literature. By now, 
static hairy black hole solutions have been found numerically for a variety of matter sources 
coupled to general relativity. Their horizon attributes, such as the dependence of the surface 
gravity on the horizon area, have been plotted and their stability has been studied in some 
detail. The literature in the field has grown steadily over the last decade with discoveries of 
new families of solutions made by combining analytic and numerical techniques. Therefore, 




(27) 




(28) 



V. DISCUSSION 
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a long list of facts regarding hairy black holes is now available. In other branches of physics, 
when sufficient data accumulates either through experiments or numerical simulations, the 
subject is generally considered to be ripe for phenomenology; simple heuristic models are 
used to put 'order' in the data. In a similar spirit, we introduced a heuristic model of hairy 
black holes as bound states of ordinary black holes and solitons and used it to account for 
qualitative features of a significant fraction of the accumulated data. 

The notion of the horizon mass played a key role both in motivating the model and in its 
applications. In particular, this notion could be used effectively to gain qualitative insight in 
to 'why' hairy black holes are unstable, and more strikingly, 'why' the magnetically charged 
Reissner-Nordstrom black holes are unstable in the Einstein- Yang- Mills theory although they 
are stable in the Einstein-Maxwell theory. Since horizon mass plays an important role in our 
phenomenology, we also investigated certain subtleties and pointed out some of its counter- 
intuitive properties. Finally although for definiteness we restricted ourselves to the Einstein- 
Yang-Mills theory in Section |T|, our phenomenological model is valid also in theories with 
more complicated matter sources and our qualitative analysis can be repeated in this more 
general context provided one takes an appropriate care of the subtleties encountered in 
Section |^ at the cross-over points. 
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